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ABSTRACT 

We present results of a statistical study of the cosmic evolution of the mass dependent major-merger rate 
since z = 1. A stellar mass limited sample of close major-merger pairs (the CPAIR sample) was selected from 
the archive of the COSMOS survey. Pair fractions at different redshifts derived using the CPAIR sample and a 
local K-band selected pair sample show no significant variations with stellar mass. The pair fraction exhibits 
moderately strong cosmic evolution, with the best-fitting function of f pa i r = 10~ L88(±0 03) (1 + z ) 2 - 2 (± - 2 ). The best- 
fitting function for the merger rate is R mg (Gyr 1 ) = 0.053 x (M star /lO ia7 M ) a3 (l +z) 2 2 /(l +z/8). This rate 
implies that galaxies of M sta r ~ 10 10_1L5 M have undergone ~ 0.5 - 1.5 major-mergers since z=l. Our results 
show that, for massive galaxies (M star > 10 10 5 M ) at z < 1, major mergers involving star forming galaxies (i.e. 
wet and mixed mergers) can account for the formation of both ellipticals and red quiescent galaxies (RQGs). 
On the other hand, major mergers cannot be responsible for the formation of most low mass ellipticals and 
RQGs of M star 10 10 - 3 M . Our quantitative estimates indicate that major mergers have significant impact on 

the stellar mass assembly of the most massive galaxies (M star > 10 1L3 M ), but for less massive galaxies the 
stellar mass assembly is dominated by the star formation. Comparison with the mass dependent (U)LIRG rates 
suggests that the frequency of major-merger events is comparable to or higher than that of (U)LIRGs. 
Subject headings: galaxies: interactions — galaxies: evolution — galaxies: starburst — galaxies: general 



1. INTRODUCTION 

Galaxy mergers have been fascinating astronomers for a 
long time, ever since they were recognized (see the review 
of Schweizer (1996)). Major mergers of galaxies of nearly 
equal mass stand out because of the more spectacular tidal 
and dynamical effects (Toomre 1978), and many nearby ma- 
jor mergers have been extensively studied (Toomre & Toomre 
1972; Whitmore et al. 1995; Hibbard & van Gorkom 1996; 
Hibbard&Yun 1999; Xu et al. 2000; Wang et al. 2004). It 
has been well documented that major mergers can induce en- 
hanced star formation (Kennicutt et al. 1987; Xu & Sulentic 
1991), trigger extreme starbursts and active galactic nuclear 
(AGN) activities (Sanders et al. 1988; Sanders & Mirabel 
1996; Dasyraetal. 2006), and transform spiral galaxies 
into elliptical galaxies (Toomre 1978; Schweizer 1982; 
Genzel et al. 2001). They dominate among the extreme 
starbursts such as luminous infrared galaxies (LIRGs, with 
SFR £ 20 M yr -1 ) and ultra-luminous infrared galaxies 

(ULIRGs, with SFR > 200 M yr 1 ; Sanders & Mirabel 

1996). On the other hand, statistically, major mergers play 
minor roles in the processes such as star formation and mass 
growth of z ~ galaxies in general. Only ~ 1 - 2 percent of 
galaxies are involved in close major mergers (Xu et al. 2004; 
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Patton & Atfield 2008; Domingue et al. 2009), and only - 2 
- 3 percents of star formation rate density (SFRD) in the z=0 
universe is due to close major mergers (Xu et al. 2010). 

Are mergers more important in the earlier universe? In- 
deed, in the hierarchical structure formation paradigm of 
the contemporary cosmology, galaxy and dark matter halo 
(DMH) merging is one of the most significant processes 
affecting the evolution of structures in the early universe, 
and is largely responsible for the growth of massive dark 
matter halos and the buildup of galaxies (Kauffmann et al. 
1993; Lacey & Cole 1993; Khochfar & Burkert 2005). Many 
observations of intermediate/high redshift peculiar galaxies 
and galaxy pairs have found strong evolution in the merger 
rate, up to (1 +z)^~6 (Brinchmann et al. 1998; LeFevre et al. 
2000; Conselice et al. 2003; Bell et al. 2006; Conselice 2006; 
Conselice et al. 2009; Kampczyk et al. 2007; Kartaltepe et al. 
2007; Rawat et al. 2008), and show evidence for mergers to 
dominate the total star formation rate in the universe of z ^ 1 

(Zheng et al. 2004; Hammer et al. 2005; Bridge et al. 2007). 
On the other hand, weak merger rate evolution (~ (1 +z) 0-5 ) 
has been found by other studies of intermediate/high red- 
shift mergers (Carlberg et al. 2000; Lin et al. 2004; Lotz et al. 
2008; Robainaetal. 2010; Manetal. 2011), and many au- 
thors have argued that at z £ 1 the SFRD is still predomi- 
nantly contributed by isolated late type galaxies rather than by 
mergers (Flores et al. 1999; Bell et al. 2005; Lotz et al. 2008; 
Jogee et al. 2009). Furthermore, mergers may no longer be 
the predominant population among LIRGs and even ULIRGs 
at z > l(Melbourne et al. 2005; Daddi et al. 2010). 

In this paper, we address the controversy on the merger 
rate evolution using data in the COSMOS field (Scoville et al. 
2007). Major reasons for the controversy include different 
definitions of major mergers, errors in the merger time scale, 
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and biases in the merger sample selections (see Xu et al. 2010; 
Hopkins et al. 2010a). Studies using merger samples selected 
from peculiar galaxies (LeFevre et al. 2000; Conselice et al. 
2003; Conselice 2006; Jogee et al. 2009; Conselice et al. 
2009) are vulnerable to contaminations due to minor merg- 
ers and to uncertainties of the time scale for the detectable 
peculiarity such as tidal tails, bridges, plumes, and other dis- 
tortions (Lotz et al. 2010). In contrast, we selected our merger 
samples from close (5 < r pro j < 20 fr 1 kpc) major-merger 
pairs (stellar mass ratios < 2.5). These pairs have reasonably 
well understood merger time scales (Kitzbichler & White 
2008; Lotz et al. 2010). 

Our pair sample is selected from the photo-z catalog of 
the COSMOS field (Ilbert et al. 2009) and is stellar mass 
limited, including massive galaxies (M star > 1O 9 M ) in the 
photo-z range of 0.2 < z p h t < 1. The COSMOS sample has 
the best photo-z 's, measured using data of ~ 30 photomet- 
ric bands covering the entire UV - infrared range, for more 
than 100,000 galaxies with nearly 100 percent completeness 
(Ilbert et al. 2009). This enables us to obtain a pair sample 
that is ~ 70% complete. By comparison the pair samples 
in the studies of Patton et al. (2002), Lin et al. (2004), and 
Bundy et al. (2009), using pairs selected from spectroscopic 
surveys, are only ~ 10 - 20% complete. Given the rather 
complex spectroscopic selection functions in those studies, 
the corrections for the incompleteness may lead to substan- 
tial uncertainties in the results. 

Photo-z selected pairs of 0.2 < z pho t < 1.2 in the COSMOS 
field were studied by Kartaltepe et al. (2007, hereafter K07). 
Their sample is different from ours in two major respects: 
(1) It is not confined to major mergers (i.e. no constraint on 
the mass ratios or luminosity ratios between primaries and 
secondaries); (2) it is a absolute magnitude limited sample 
(My < -19.8). With more rigorously defined major-merger 
pair samples, our study shall improve upon the results of K07. 
Also, with well determined stellar mass for every galaxy in 
the sample, we shall study the mass dependence of the merger 
rate evolution. 

We will confine our analysis to galaxy pairs of z < 1, be- 
cause the photo-z 's and stellar mass estimates of z > 1 galax- 
ies are less accurate (Ilbert etal. 2010). We will derive the 
merger rates for galaxies of different stellar masses and red- 
shifts, and separate dry mergers (E+E pairs) and wet/mixed 
mergers (S+S and S+E mergers). The selection of COSMOS 
pairs in the redshift range of 0.2 < z < 1 is described in Sec- 
tion 2. Corrections for the incompleteness and for contami- 
nations by spurious pairs are presented in Section 3. A z=0 
pair sample, which sets the local benchmark for the evolution 
study, is presented in Section 4. The mass dependent merger 
rate and its evolution since z=l are presented in Section 5. 
The contributions of major mergers to the galaxy assembly 
and elliptical galaxy formation since z=l is investigated in 
Section 6, and comparisons to the mass and redshift depen- 
dence of the (U)LIRGs' abundance is in Section 7. Section 8 
is devoted to a summary of our main results. Throughout this 
paper, we adopt the A-cosmology with Q m = 0.3 and ft \ = 0.1, 
and H = 70 (km sec _1 Mpc _1 ). 

2. THE COSMOS PAIR (CPAIR) SAMPLE 

We selected major-merger pair candidates using a parent 
sample of galaxies constructed from that used by Drory et al. 
(2009, hereafter D09) in their study of galaxy stellar mass 
function (GSMF), which is in turn selected from the COS- 
MOS photo-z catalog (Ilbert et al. 2009) using the following 



Table 1 

Parent Sample 



z min 


z max 


Volume 


log(M min ) 


Number of Galaxies"^ 






(10 6 Mpc 3 ) 


(M ) 


SFGs 


RQGs Total 


0.2 


0.4 


0.56 


9.0 


6787 


2039 8826 


0.4 


0.6 


1.23 


9.4 


6169 


1526 7695 


0.6 


0.8 


1.92 


9.8 


6745 


1981 8726 


0.8 


1.0 


2.53 


10.2 


6610 


2287 8897 





Note. 



: f Number of galaxies with log(M sta r) > log(M m i n ). 

criteria: 0.2 < z < 1, K s < 24 and /^b < 25 - 5 - 

The D09 sample has 138001 galaxies, divided into four 
photo-z bins of width of Az = 0.2, and into "active" (star 
forming galaxies, SFGs) and "passive" (red quiescent galax- 
ies, RQGs) populations according to the SED type of the best 
fitting template (Ilbert et al. 2009). The stellar mass of galax- 
ies, M star , is derived through a stellar population synthesis 
model fitting (the Chabrier IMF), using the photo-z and pho- 
tometric data in the u* (CFHT), Bj, V 7 , g + , r + , i + , z + (Sub- 
aru), / (UKIRT) and K s (CFHT) bands. Typical uncertainties 
in M star is between 0.1 dex and 0.3 dex at 68% confidence 
level, depending on spectral type and the S/N of the photom- 
etry (D09). In the four photo-z bins, the completeness limits 
for SFGs and RQGs are log(M star /M ) = [8.3,8.9,9.2,9.4] 
and log(M star /M ) = [8.9,9.2,9.8, 10.1], respectively. 

For galaxies in our parent sample, the stellar mass is taken 
from D09. In order to be complete for both SFGs and RQGs, 
we imposed a mass limit (log(M m i n )) on each of the photo-z 
bins in the parent sample. The mass limit and the number of 
galaxies above the limit are listed in Table 1. There are 34144 
galaxies in the parent sample. 

The pair sample is also divided into four photo-z bins. The 
selection criteria are: 

: (1) The primary galaxy has log(M star ) > log(Mii m ), with 
log(Mi im /M ) = [9.4,9.8, 10.2, 10.6] for the four red- 
shift bins, respectively. Mii m 's are 0.4 dex above the 
M m i n 's of the parent sample (Table 1). 

: (2) the difference in M star between the primary galaxy and the 
secondary galaxy is less than 0.4 dex: Alog(M star ) < 
0.4. 

: (3) the redshift difference between the two components, 
Azphot = |z|L- z pEotl> satisfies Az phot /(l+zp < 0.03. 

: (4) the projected physical separation (r pro j) is in the range of 
5 < r proj < 20 hr 1 kpc. 

Compared to the selection criteria for local pairs described 
in Xu et al. (2004), we replaced the rest- frame K-band se- 
lection by a stellar mass selection in criteria (1) and (2). 
The Spizer-IRAC 3.6/im and 4.5 jim bands, which encom- 
pass the rest- frame K-band emission for galaxies of 0.6 ^ z < 

1, have relatively low angular resolution compared to the 
HST/ground-based optical and NIR data. Using IRAC data 
would have resulted in larger confusion errors in the stellar 
mass for pairs with separation £ 2" . At the same time, it 
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Table 2 

Sample of Paired Galaxies in COSMOS 



z min 


z max 


Number of Galaxies ^ 




log(Mi im ) 




Number of Galaxies with M 


5 tar > Mum* 








in iso. 


in multi. 


Total 


(M ) 


in iso. 


in multi. 


Total 


SFGs* 


RQGs* 






pairs 


systems 






pairs 


systems 








0.2 


0.4 


144 


6 


150 


9.4 


128 


5 


133 


78 


55 


0.4 


0.6 


100 


3 


103 


9.8 


93 


3 


96 


61 


35 


0.6 


0.8 


144 


22 


166 


10.2 


126 


20 


146 


109 


37 


0.8 


1.0 


174 


24 


198 


10.6 


131 


21 


152 


81 


71 





Note. 

: f Including all primaries and all secondaries. 

: i Including all primaries and those secondaries with log(M star ) > log(Mii m ). 

: * SFGs ("active galaxies") and RQGs ("passive galaxies") classifications were taken from D09. 



was shown in D09 that for field galaxies of z < 1 the stellar 
mass derived using the HST/ground-based optical and NIR 
data are nearly identical to those derived using data including 
the IRAC fluxes (Ilbert et al. 2010). 

Criteria (1) and (2) guarantee that our pair sample is not 
affected by the "missing secondary" bias (Xu et al. 2010). 
Criterion (3) is set to minimize the incompleteness due to the 
photo-z error while ensure that the reliability of the sample is 
not significantly compromised. This issue will be addressed 
in detail in Section 3. 

Using these criteria, 417 pair candidates were selected. In 
order to exclude spurious pairs due to imaging artifacts, visual 
inspections were carried out on the HST-ACS images (F814 
band, Koekemoer et al. 2007). Among 417 pair candidates, 
335 were covered by the HST survey. All pair candidates out- 
side the area of HST-ACS imaging (82/417 = 19.7%) were 
dropped from the final pair sample. Fourteen spurious pairs 
were identified: 8 have wrong astrometry (i.e. no source ap- 
pears at the sky coordinates) for at least one of the two galax- 
ies and 6 are pieces of single large disc galaxies. 

Some galaxies were found repetitively in multiple pair can- 
didates. These 40 pair candidates consist of 17 triplets and 
1 quartet, including 55 galaxies. The remaining 281 can- 
didates are isolated pairs, including 562 galaxies. Because 
galaxies in triplets/quartets are less than 10% of the total sam- 
ple, we will not distinguish them from paired galaxies. Our 
final COSMOS Pair sample (hereafter CPAIR) includes 617 
paired galaxies, found in both isolated pairs and multiple sys- 
tems. Among them, 527 (including both primaries and secon- 
daries) have log(M star ) > log(Mii m ), and the remaining (90, all 
secondaries) have (log(Mi im )-0.4) < log(M star ) < log(Mi im ). 
Statistics of the sample are listed in Table 2. 

3. CPAIR SAMPLE: INCOMPLETENESS AND SPURIOUS PAIRS 
FRACTION 

Much of the discrepancies between different results on 
merger rate evolution can be attributed to various biases caus- 
ing incompleteness (missing of true mergers) and contamina- 
tions of spurious mergers in merger samples. Therefore it is 
important to investigate thoroughly all such biases and correct 
them in merger statistics. 

3.1. Incompleteness due to Missing Very Close Pairs 

Photometric data of the photo-z catalog (Ilbert et al. 
2009) were obtained using the SExtractor in dual mode 



(Bertin & Arnouts 1996). Images in all bands were degraded 
to a common PSF of FWHM = 1 .5", and the photometry was 
done with a constant aperture of r = 1.5" (Capak et al. 2007). 
Because of the limited angular resolution of the photo-z cat- 
alog, very close pairs with angular sepatation e £ 2" are in- 
complete in the pair sample. Exploiting the COSMOS HST- 
ACS lensing catalog (Leauthaud et al. 2007, 2010), we esti- 
mated this incompleteness to be [0.01, 0.06, 0.08, 0.20] for 
the four redshift bins, with no significant mass dependence. 
The full analysis can be found in Appendix A. 

3.2. Incompleteness due to Photo-z Errors and Spurious 

Pairs due to Projection 

For pairs of e > 2.0", major cause of the incompleteness is 
due to photo-z errors, which have a non-negligible probabil- 
ity of being so large that a real pair with Av < 500 km sec -1 
(corresponding to Az/(l+z) < 0.0017) can have a mea- 
sured Az p hot/(l + z p hot) > 0-03 and therefore be missed by 
the CPAIR sample. Also, the photo-z selection criterion and 
the photo-z errors can introduce spurious pairs whose velocity 
difference Av is larger than 500 km sec -1 . Using Monte Carlo 
simulations, we estimated the incompleteness and the spuri- 
ous pair fraction (hereafter SPF) to be [0.21, 0.21, 0.23, 0.25] 
and [0.07, 0.08, 0.10, 0.09], repectively, for the four refshift 
bins. The full analysis is presented in Appendix B. 

3.3. Clustering Effect on Spurious Pair Contaminations 

Bell et al. (2006) found that in the COMBO-17 survey, the 
projected two-point correlation functions of massive galax- 
ies with 0.4 < Zphot < 0.8 can be well fitted by a power-law 
w(r) oc r~ 7 down to r = 15 kpc, with the value of the power- 
index 7 consistent with 2. Based on this result (see also 
Robainaetal. 2010), we made a simple estimation for the 
effect of galaxy clustering that was neglected in the Monte 
Carlo simulations. It should be pointed out that, different from 
Patton et al. (2000) and Bell et al. (2006), we assumed that 
the boundary separating physical (i.e. gravitationally bound) 
and unphysical pairs, r p , is much larger than ri = 20 h _1 kpc, 
the separation limit in our pair selection. This is because the 
merger time scales derived by Lotz et al. (2010), as adopted 
in this work (Section 5.2), are for close pairs with projected 
separation r < 20 h -1 kpc. 

For a given two-point correlation function £(r) = (ro/r) 7 , the 
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Table 3 

Completeness and Reliability Corrections for CPAIR Sample 



z min 


z max 




Completeness Correction 






Reliability Correction 








due to mis sin 


g due to photo-z 


Combined 


due to random 


due to 


due to pairs ot 


Combined 






very close pairs errors 




projection 


clustering 


Av > 500km/sec 




0.2 


0.4 


0.99 ±0.01 


0.79 ±0.05 


0.78 ±0.05 


0.93 ±0.03 


0.94 ±0.05 


0.91 ±0.03 


0.80 ±0.06 


0.4 


0.6 


0.94 ±0.03 


0.79 ±0.05 


0.74 ±0.05 


0.92 ±0.03 


0.94 ±0.05 


0.91 ±0.03 


0.79 ±0.06 


0.6 


0.8 


0.92 ±0.03 


0.77 ±0.05 


0.71 ±0.05 


0.90 ±0.03 


0.94 ±0.05 


0.91 ±0.03 


0.77 ±0.06 


0.8 


1.0 


0.80 ±0.04 


0.75 ±0.05 


0.60 ±0.05 


0.91 ±0.03 


0.94 ±0.05 


0.91 ±0.03 


0.78 ±0.06 



Figure 1. Illustraction for Eq. 1. ri is the maximum projected separa- 
tion in the pair selection, r p the outer boundary of physical (i.e. gravi- 
tationally bound) pairs, and ta the distance range along the line-of-sight 
that corresponds to the pair selection criterion for the photo-z difference: 
| Azphot|/(l + z phot) < 0.03. The shadowed areas are places where spurious 
companions are located. 

additional SPF (77) due to clustering can be estimated as fol- 
lows: 



7]-- 



n(n/r) 



sin(0)dO 



4irn 



rarcsin{r\ jr) 



'sin(6)d9 
(1)' 

The numerator on the right-hand side of the equation is the 
probability to find spurious companions near a galaxy in both 
foreground and background (in the shadowed areas in Fig. 1). 
Here ta is ~ 100 Mpc, corresponding to the pair selection 
criterion of | Az p h t|/(1 +Zphot) < 0.03. The denominator is 
the probability of finding both real and spurious companions 
with a projected separation of r < ri. The relations between 
parameters ri, r p and ta are illustrated in Fig. 1. For 7 = 2 and 
ri 



<< r p << ta, Eq. 1 can be approximated by: 



77 = 



n/r p 



(2) 



Assuming r p = 100 fr 1 kpc, Eq. 2 gives 77 = 0.06. It is worth 
noting that: (1) 77 is comparable to the SPF found by the 
Monte Carlo simulations for the random associations; (2) 77 
is constant against the redshift; (3) given the uncertainties in 
r p and 7, we shall assume an relatively large error of 0.05 for 
77; (4) this correction also applies to pairs selected spectro- 
scopically because the condition r a > > r p is still valid even 
for the selection criterion of Av < 500 km sec -1 . 

3.4. Fraction of Physical Pairs with Av > 500 km sec -1 

We excluded physical pairs with 500 < 
Av ;$ 1000 km sec -1 from the merger rate analysis since 

they have very uncertain and long merger time scales. These 
high Av pairs, found in the group/cluster environments 
(Domingue et al. 2009), are not included in the above 
estimates for projected unphysical pairs. Here we make 
a separate correction for them. Using local pairs taken 
from Domingue et al. (2009), selected using nearly identical 
selection criteria as the CPAIR sample (cf. Section 4) except 
for that spec-z were used and they include all pairs (isolated 



or in groups/clusters) with Av spec < 1000 km sec \ we find 
that the fraction of physical pairs with Av > 500 km sec -1 is 
9.6 ±3.0% (Appendix C). 

3.5. Completeness and Reliability of the CPAIR Sample 

In Table 3 we listed estimates of the correction factors for 
subsamples in the four redshift bins. The combined complete- 
ness correction factor is the product of that due to missing 
very close pairs and that due to photo-z errors. The combined 
reliability correction factor, defined as 1 - SPF, is the product 
of that due to random projections, the additional correction 
due to the clustering effect, and that due to the contamination 
of physical pairs with Av > 500 km sec -1 . The combined 
completeness correction factor varies in the range of 0.60 - 
0.78 between the 4 redshift bins. The combined reliability 
correction factor (~ 0.79) is rather constant against the red- 
shift. 

As an independent check, exploiting a sample of spec-z 
pairs (the ZPAIR sample) selected from the zCOSMOS sur- 
vey (Lilly et al. 2007), we made an empirical analysis on the 
completeness and reliability of the CPAIR sample. This re- 
sulted in an estimate of 0.86 ± 0. 12 for the completeness cor- 
rection factor due to photo-z errors (the incompleteness due to 
missing very close pairs cannot be checked with spec-z pairs), 
which is consistent (within 1-cr) with the result of the Monte 
Carlo simulations (Table 3). The estimate for the reliability 
correction is also 0.86 ±0.12, again consistent (within l-a) 
with the values of the combined reliability correction in Ta- 
ble 3. The details of the analysis are presented in Appendix D. 

4. LOCAL PAIR SAMPLE 

Pair statistics in the local universe were carried out using 
an updated version of the KPAIR sample by Domingue et al. 
(2009), a close major-merger pair sample selected in the Ri- 
band from cross matches between 2MASS and SDSS-DR5 
galaxies. The update includes following modifications: 

: (1) Stellar masses of galaxies in KPAIR and in its parent 
sample are multiplied by a factor of 10~ 39 . This is 
the average difference between the mass estimated us- 
ing the total K s band luminosity and a Salpeter IMF 
(Domingue et al. 2009), and the mass estimated using a 
Kroupa IMF by Kauffmann et al. (2003). Because the 
mass estimated using the Kroupa IMF and that using 
the Chabrier IMF are nearly identical (Kauffmann et al. 
2003), this modification makes the masses in the local 
KPAIR sample and those in the CPAIR sample consis- 
tent. 

: (2) In order to avoid possible bias due to the local over- 
density (associated with the local super-cluster), a 
lower redshift cut-off of v > 2000 km sec -1 (z > 
0.0067) is introduced. 
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Table 4 

Mass Dependent Pair Fractions 



Mass Bin 


z = 




0.2 < z < 0.4 


0.4 < z < 0.6 


0.6 < z 


< 0.8 


0.8 < z < 1 




N pg /N G 


fpair 


N Pg /N G 


fpair 


N Pg /N G 




N Pg /N G 




N Pg /N G 


9.4 < log(M) < 9.8 


0.014±0.011 


21 126 


0.027 ±0.006 


44/2034 




0/0 




0/0 




0/0 


9.8 < log(M) < 10.2 


0.011 ±0.006 


7/524 


0.020 ± 0.006 


24/1482 


0.030 ± 0.006 


44/1977 




0/0 




0/0 


10.2 < log{M) < 10.6 


0.013 ±0.002 


41/2775 


0.043 ±0.010 


39/1174 


0.030 ±0.007 


34/1524 


0.038 ±0.007 


71/2542 




0/0 


10.6 < log(M) < 11.0 


0.014 ±0.002 


97/5826 


0.035 ±0.011 


19/ 706 


0.025 ±0.007 


18/ 955 


0.044 ±0.008 


62/1913 


0.045 ±0.010 


110/3229 


11.0 < log(M) < 11.4 


0.011 ±0.002 


68/4520 


0.042 ±0.020 


6/183 


(0) 


0/267 


0.021 ±0.008 


10/ 645 


0.061 ±0.015 


39/1051 


11.4 < log(M) < 11.8 


0.014 ±0.007 


8/417 


(0) 


0/17 


(0) 


0/29 


(0) 


0/65 


0.048 ±0.029 


3/ 103 


total 


0.013 ±0.001 


223/14188 


0.030 ±0.010 


132/5596 


0.026 ±0.007 


96/4752 


0.038 ±0.009 


143/5165 


0.056 ±0.011 


152/4383 



(3) Pairs with 500 < Av < 1000 km sec -1 are excluded. Most 
of these pairs are in cluster environment, and may not 
be gravitationally bound. Excluding them improves the 
accuracy of the merger rate estimate. 



(4) The magnitude limit is set at K s 
limit of the KPAIR sample. 



12.5, the completeness 



There are 18,081 galaxies in the parent sample that are 
brighter than K s = 12.5, of which 14,813 have measured red- 
shifts (redshift completeness of B z _ com p = 0.82), and 14,218 
are in the range of 0.0067 < z < 0.1. The new paired galax- 
ies sample has 221 galaxies, all brighter than K s = 12.5. 
Among them 188 are in pairs with two measured redshifts 
of Av spec < 500 km sec -1 . These redshifts are in the range of 
0.0067 < z < 0.1, with the median of z=0.042. The remain- 
ing 33 galaxies are in single-redshift pairs (i.e. only one of 
the component galaxies having measured redshift). 

5. MASS DEPENDENT MERGER RATES 

5.1. Pair fraction 

The local pair fraction is calculated using the following for- 
mula: 



f, 



pair,0 — 



Bz-comp X (1 -7?) ^ N 2z +Ni z X (1 - Qsp urious) 

1 w/ 



A 



N G ,o 



where B 



z-comp 



= 0.82 is the redshift completeness of the par- 
ent sample, Ao = 0.89 the completeness of the local pair sam- 
ple (Domingue et al. 2009), (1 -rj) = 0.94 the clustering re- 
lated reliability correction factor found in Section 3.3. N 2z 
and Niz are numbers of galaxies in pairs of two measured 
redshifts and single redshifts, respectively, Nq,o the num- 
ber of galaxies in the parent sample with measured redshift 
(in the range of 0.0067 < z < 0.1). Q spu rious = 0.2 is the 
probability for a single-redshift pair to be a spurious pair 
(Domingue et al. 2009). Note that Ao ^Bz_ comp because pairs 
with single measured redshifts were included. Also, we sup- 
plemented the SDSS redshifts of paired galaxies with red- 
shifts found in the literature and in our own redshift obser- 
vations (Domingue et al. 2009). 

The pair fractions in the COSMOS field are estimated as 
follows: 

f . - Q reli x ^Ei (A) 

pair "C comp x(l-D ACS ) N G ' 

where Q reli = [0.80,0.79,0.77,0.78] is the reliability and 
Ccomp = [0.78,0.74,0.71,0.60] the completeness (for the four 
photo-z bins) of the pair sample, respectively; D A cs = 0.197 
is the fraction of pair candidates without ACS images (not 
included in the final pair sample, see Section 2); N pg is the 



number of interacting galaxies in the pair sample, and Nq the 
number of galaxies in the parent sample. 

The relative error of the pair fraction (i.e. error/f pa i r ) can be 
estimated as the quadratic sum of the random error cr rms and 
the cosmic variance cr var i 7 : 



2 , 2 



The random (binomial statistics) error is 



1-f 



pair 



N 



(5) 



(6) 



Pg 



For the local sample, we adopted the approximation N pg = 
N2z + Niz(l - Qspurious)- The cosmic variance is given by 
(Peebles 1980; Somerville et al. 2004): 



cr vari =J 2(7)x(roA: 



samp 



(7) 



where ro and 7 are the parameters in the two-point correla- 
tion function £(r) = (ro/r) 7 , r samp the radius of the sampling 
volume, and J2 a function of 7: 



h = 



72 



(3-7)(4- 7 )(6-7)2^ 



(8) 



The correlation function parameters for local galaxies of dif- 
ferent masses were taken from Zehavi et al. (2005). For 
galaxies of z > 0.2, we assumed 7 = 1.8 and derived the 
ro values by interpolating the measurements for galaxies 
of different masses and redshifts by Zehavi et al. (2005), 
Meneux et al. (2008) and Foucaud et al. (2010). Cosmic vari- 
ances for the integral pair fractions (i.e. not divided into mass 
bins) were taken from that for galaxies in the mass bin of 
10.6<log(M star /M o )<ll. 

The pair fractions so calculated are listed in Table 4 and 
plotted in Fig. 2. For local galaxies, we confirm the con- 
clusion of Domingue et al. (2009) that there is no significant 
mass dependence of the pair fraction. The integral pair frac- 
tion at z=0 is 1.3 ±0.1%, slightly lower than the result of 
Domingue et al. (2009) which is 1.6 ±0.1%. The difference 
is due to two factors: (1) the sample of Domingue et al. (2009) 
includes all pairs of Av < 1000 km sec -1 whereas local pairs 
in this work are restricted to pairs of Av < 500 km sec -1 ; (2) 
the application of the clustering related reliability correction 
factor ((I-77) = 0.94). 

There is no evidence for a significant mass dependence of 
the pair fractions in higher photo-z bins, either. The trend 

7 The pair fraction, f pa j r , is proportional to the probability of finding a 
second galaxy within a spatial separation r from a given galaxy: P(r) = 
47rn J* [I + £(r)]r 2 dr, where n is the number density of galaxies and £ the 
two-point correlation function. Hence f pa i r is proportional to n, and therefore 
the cosmic variance has the same effect on f pa i r as on the number density. 
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Figure 2. Mass-dependent pair fractions in photo-z bins. 



0.6 
0.5 

^ 0.4 

S 0.3 



0.2 



0.1 



0.0 



$ 0.4<z<0.6 ij 0.8<z<1.0 
$ 0.2<z<0.4 4 0.6<z<0.8 
f z=0 




9.0 9.5 10.0 10.5 11.0 11.5 12.0 

log (M sta /M sun ) 

Figure 4. Mass-dependent differential merger rate in different photo-z bins. 
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Figure 3. Pair fraction evolution. The solid line is the least-square fit to 
total pair fraction (pair fraction of all galaxies regardless of the stellar mass) 
vs. redshift relation, specified by f pair = 10- L88 (± 003 )(1 + z ) 2 - 2 (± - 2 ). The 
dashed line is the least-square fit to the pair fractions in the mass bin of 10.6 < 
log(M st ar/M ) < 11, specified by f pair = io- L84 (±°- 05 )(l + z ) 2 - 1(±a3) . 



for pair fractions to increase with redshift can be seen in all 
mass bins, though with substantial scatter. The major reason 
for the large scatter is the cosmic variance, given the rela- 
tively small volume explored by the COSMOS survey in each 
photo-z bin (Table 1). In particular, there is a strong den- 
sity enhancement in the photo-z bin of 0.2 < z < 0.4 (D09), 
which biases the pair fraction towards a higher value (f pa i r is 
proportional to the density). Cosmic variance often dominates 
the total error in the pair fraction: its contribution is usually 
> 50% except in those bins where the number of paired galax- 
ies, N pg , is less than 10 (hence the random error is large). In 
Fig. 3, we plot the redshift dependence of the integral pair 
fraction (pair fraction of all galaxies regardless of the stellar 
mass), and that of the pair fraction for galaxies in the mass bin 
10.6 < log(M sta r/M ) < 1 1 (the bin encompassing the M* tar ). 
These two results are very close to each other, in agreement 
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Figure 5. Comparisons of mass-dependent pair fractions of this work with 
those of Bundy et al. (2009). Both results are converted to the f pa i r for merg- 
ers of mass ratio, /x, < 3 (see text). 

with our conclusion that the pair fraction does not vary signif- 
icantly with stellar mass. 

The least- square fit to the redshift dependence of the inte- 
gral pair fractions is f pair = lO" 1 - 88 ^ - 03 ^ +z) 2 - 2(±a2) . For the 
pair fractions in the mass bin of 10.6 < log(M star /M ) < 11, 
the best fit is f pair = lO- 1 - 84 ^ - 05 ^ + z) 2 - 1 < ±0 - 3 >. 

5.2. Differential Major-Merger Rate 

The differential major-merger rate is the probability for 
each galaxy to be involved in a major merger per Gyr: R mg oc 
fpair /T mg , where T mg is the merger time scale in Gyr. Because 
the physical process of a galaxy merger is very complex (see 
Hopkins 2010b for a review), T mg has been a major source 
of uncertainty in the merger rate studies. In the literature the 
most common approach has been the approximation of T mg by 
the dynamical friction time scale (Binney & Tremaine 1987; 
Pattonetal. 2000; Jiang etal. 2008; Kitzbichler & White 
2008). Kitzbichler & White (2008) studied the T mg in a semi- 
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Table 5 

Mass Dependent Differential Merger Rate 



Mass Bin 






RmgCGyr" 1 ) 








z = 


0.2 < z < 0.4 


0.4 < z < 0.6 


0.6 < z < 0.8 


0.8 < z < 1 


9.4 < log(M) < 9.8 
9.8 < log{M) < 10.2 
10.2 < log(M) < 10.6 
10.6 < log(M) < 11.0 
11.0 < log(M) < 11.4 
11.4 < log(M) < 11.8 


0.023 ±0.019 
0.025 ±0.012 
0.037 ± 0.008 
0.055 ± 0.007 
0.062 ±0.011 
0.096 ±0.046 


0.050 ±0.011 
0.049 ±0.013 
0.132 ±0.030 
0.140 ±0.042 
0.226 ±0.107 


0.068 ±0.014 
0.089 ± 0.020 
0.099 ± 0.029 


0.1 12 ±0.020 
0.171 ±0.033 
0.108 ±0.040 


0.211 ±0.037 
0.302 ±0.076 
0.313 ±0.191 





analytical model built on the results of the Millennium Simu- 
lation. They assumed circular orbits to estimate the dynamical 
friction process, and found relatively weak mass and redshift 
dependence in the form of T mg oc M~^ r 3 (l +z/8). However, 
their dynamical friction time is only appropriate for small 
satellite galaxies at large radii. For the massive close major- 
merger pairs in our samples, it becomes a very poor approxi- 
mation because of two issues (Hopkins et al. 2010a): (1) An- 
gular momentum loss at these radii is not dominated by the 
dynamical friction, but rather by exchange in strong reso- 
nances between the baryonic components that act much more 
efficiently. (2) By these radii, even the initially circular orbits 
have become highly radial, leading to shorter merger times. 

On these reasons, we instead estimated T mg using the re- 
sults of Lotz et al. (2010). They carried out high resolution 
hydro-dynamical simulations for a large number of galaxy 
mergers with diverse initial conditions, and derived T mg at 
different projected separations in a view-angle averaged for- 
mat. These simulated mergers have line-of- sight velocity dif- 
ference Av < 500 km sec -1 , identical to the pairs in our sam- 
ple. Nine mergers in Lotz et al. (2010) have baryonic mass 
ratios < 3 (i.e. one 1:1 merger and eight 3:1 mergers) with 
M st ar > 10 10 - 7 M . In the bin of 5 < r pro j < 20 IT 1 kpc, 
these nine mergers have an average merging time scale of 
T mg = 0.30 ± 0.06 Gyr. The three 1:1 mergers in Lotz et al. 
(2010) have stellar masses in the range of 10 9 7 - 10 10 7 M , 
and their T m g s in the 5 < rproj ^ 20 h * kpc bin show a 
weak mass dependence of T mg ~ M~ t ^ r 2 , consistent with what 
found by Kitzbichler & White (2008). Because Lotz et al. 
(2010) did not study the red- shift dependence of T mg , and 
because their sample is too small (3 mergers) to derive a 
meaningful mass dependence, we adopt the relation found 
by Kitzbichler & White (2008), namely T mg ex M-° ar 3 (l +z/8). 
These dependences are sufficiently weak that the associated 
uncertainties will not have any significant effect on our re- 
sults. The final merger time scale we adopted is: 



T mg = 0.3 Gyr x 



1O 1O - 7 M 



-0.3 

('♦!)■ 



(9) 



And the differential merger rate for major mergers of mass 
ratio < 3 is: 

R mg = A X fpair/T m g7 (10) 

where A = 1 . 19 is the factor converting the pair fraction in this 
work (for mergers of mass ratio /i < 2.5) to that of mergers 
of mass ratio /i < 3. Here we assumed that f pa i r oc log(/i max ) 
(Appendix E). In Table 5 and Fig. 4 we present our results on 

Rmg- 

Using the least-square fit to the pair fraction evolution, 
f pair = i - 1 - 88 (± 03 )(l +z )2-2(±o.2) 9 and Eq 9? we derived a best _ 



fit function for the mass-dependent R mg evolution: 



0.3 



R mg (M star ,z) = 0.053 x J ^ (Gyr-). 

(11) 

Integrating this merger rate over time, we find that the proba- 
bility for individual galaxies to be involved in a major merger 
since z=l is O.8x(M star /lO ia7 M ) - 3 . Accordingly, on aver- 
age, massive galaxies of M star ~ 1O 1O_1L5 M have undergone 
~ 0.5 - 1.5 times mergers since z=l. 

5.3. Comparisons with Previous Results 

Our results on the mass independence of the local pair 
fraction (filled circles in Fig. 2) are in agreement with those 
of Domingue et al. (2009) and Patton & Atfield (2008) while 
contradicting Xu et al. (2004), the latter were derived using 
a small sample of 19 pairs. For pairs of higher redshifts, 
Bundy et al. (2009) found a trend of positive mass depen- 
dence, which was not confirmed by our results. In Fig. 5, 
our results are compared to those of Bundy et al. (2009). In 
order to compensate the difference in the mass ratios in the 
two works (/i < 2.5 in this work and fi < 4.0 in Bundy et al. 
2009), both results are converted to the f pa i r for fi < 3 mergers. 
In Appendix E, it is shown that the pair fraction f pa i r increases 
proportionally with log(/i max ). Accordingly, the pair fractions 
from this work were scaled up by a factor of 1 . 19 = log(3) / 0.4, 
and those of Bundy et al. (2009) were scaled down by a fac- 
tor of 0.80 = log(3)/0.6. The results of Bundy et al. (2009) 
might have suffered from large uncertainties: Those obtained 
using their "method I" (projected pairs without any redshift 
information for the companions), which are plotted in Fig. 5, 
were based on pair samples with high contaminations (~60 - 
70%) of unphysical pairs; and those from their "method II" 
(spectroscopic and/or photometric redshifts for both compo- 
nents) were based on small samples (3 to 15 paired galax- 
ies in each mass/redshift bin), de Ravel et al. (2009) claimed 
evidence for strong mass dependence of the evolutionary in- 
dex of the pair fraction, in the sense that low mass pairs have 
strong pair fraction evolution (m = 3 . 1 3 ± 1 .54) and high mass 
pairs have weak evolution (m = 0.52±2.07). However their 
results are very uncertain, as indicated by their large errors, 
de Ravel et al. (2009) also claimed evidence for strong evolu- 
tion (m = 4.73 ± 2.01) in optically faint pairs (M B < -18 -Qz, 
Q=l.ll) and for weak evolution (m = 1.50 ±0.76) in opti- 
cally bright pairs (M B < -18.77 -Qz). But the low evolu- 
tionary index of the bright pairs was obtained only when they 
included in their fit the z=0 pair fraction of de Propris et al. 
(2007), one of the highest local pair fraction in the literature 
(Fig. 6). Indeed, when being calculated in the same way as 
for the evolutionary index of faint pairs (i.e. fitting only high 
z data points), the index of bright pairs is m = 3.07 ± 1.68, 
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consistent with that for faint pairs. 
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Figure 6. Comparisons of observed pair fractions in the literature. When it 
is appropriate, results of different authors were corrected so they are consis- 
tent with a common definition of close major-merger pairs with the maximum 



projected separation of r pj 



= 20 h kpc and the maximum primary-to- 



secondary mass ratio /i max = 3. Results of pair samples of different r pro j 5ir 
were corrected by assuming f pa j r oc (r pro j,max) 3_7 with 7 = 2 (Bell et al. 2006). 
For example, results of Lin et al. (2008) were divided by a factor of 1 .5 
because they had r pro j ?max = 30 h _1 kpc. The results of Bell et al. (2006) 
had r m 



= 30 kpc, corresponding to r p] 



: 21 IT 1 kpc for h = 0.7, 



very close to the r p] 



: = 20 h kpc and therefore no correction was ap- 



plied. Results of pair samples of different mass ratio limits were corrected 
by assuming f pa j r oc log(/i m ax) (Fig. E-l). These include results of this work 
(log(/imax) = 0.4), Bundy et al. (2009) and de Ravel et al. (2009) (both having 
log(/imax) = 0.6), and Patton & Atfield (2008) (log(/x max ) = 0.3). This correc- 
tion was not applied to results for pair samples without the mass ratio cutoff 
(e.g. K07, Bell et al. 2006, and de Propris et al. 2007). 

In Fig. 6, our results on the cosmic evolution of the in- 
tegral pair fraction are compared with those taken from the 
literature. It shows that the evolutionary rate derived from 
our results is in between those for the strong evolution (e.g. 
the result of K07) and for weak evolution (e.g. the result of 
Lin et al. 2008), respectively. Actually, our pair fractions in 
the photo-z bins of z = 0.2-1.0 agree well with those of K07 
in the same redshift range. The marginally significant differ- 
ence between the two evolutionary rates is mainly due to the 
relatively low pair fraction at z=0, which may be caused by 
an incompleteness associated to the "missing secondary" bias 
(Xu et al. 2004), and the relatively high pair fraction at z=1.3 
in the results of K07. We clearly see much stronger evolution 
than that of Lin et al. (2008). Their results, based on spec- 
troscopically confirmed pair samples in incomplete redshift 
surveys, may have relatively large statistical uncertainties be- 
cause of the large correction factors (a factor of > 3) for the 
incompleteness. The color based pre-selection of their red- 
shift surveys may indeed introduce biases in the pair selection, 
given the significant influence of galaxy-galaxy interaction on 
optical colors (Larson & Tinsley 1978). 

In the literature, pair fractions are often compared to 
merger fractions estimated using counts of peculiar galax- 
ies (Conselice et al. 2003, 2009; Lotz et al. 2008; Jogee et al. 
2009). In general the latter are higher than the former, be- 
cause (1) Contaminations from irregular galaxies (Jogee et al. 
2009); (2) morphologically selected merger samples based 
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Figure 7. The redshift dependence of T = 1 /R m g. The values of morpholog- 
ical based merger studies of Conselice et al. (2009) and of Lotz et al. (2008) 
were taken from Fig. 7 of Conselice et al. (2009). R m g's of Bundy et al. 
(2009) were scaled down by a factor of 0.80 = log(3)/0.6 to make it con- 
sistent with the merger rate for /1 < 3 mergers. The solid line is the inverse of 
Eq. 1 1 of this work. The dashed line is the best fit of Conselice et al. (2009). 
The dotted line is the model prediction of Hopkins et al. (2010b). 
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Figure 8. Data points with error bars are mass-dependent volume merger 
rates in different redshift bins. The dashed line is the average growth rate of 
ellipticals between z=0.3 - 0.9 and dotted line that of RQGs, estimated using 
the mass function of Ilbert et al. (2010). The solid line is the average volume 
merger rate calculated using Eq. 12 by replacing R mg with its best fit (Eq. 11) 
and averaged over the same redshift range of 0.3 - 0.9. 

on the G - M 20 method (Lotzetal. 2008, 2010) include 
minor mergers; (3) the merger time scales for morpholog- 
ically selected merger samples based on the CAS method 
are longer than the merger time scales of close major-merger 
pairs (Conselice et al. 2009). Given the different merger time 
scales for close pairs and for peculiar galaxies, it is more 
appropriate to compare the differential merger rates R mg . 
In Fig. 7 we compare the inverse of the R mg , T = 1/R m g 
(Conselice et al. 2009), of morphologically selected mergers 
by Conselice et al. (2009) and by Lotz et al. (2008) with that 
of M* ar galaxies (log(M* ar /M ) - 10.8, Ilbert et al. (2010)) 



Galaxy Pairs in COSMOS - Merger Rate Evolution Since z = 1 



Table 6 

Mass Dependent Volume Merger Rate 



Mass Bin 




log(Rv) (Mpc^dex^Gyr -1 


) 






z = 


0.2 < z < 0.4 


0.4 < z < 0.6 


0.6 < z < 0.8 


0.8 < z < 1 


9.4 < log(M) < 9.8 
9.8 < log{M) < 10.2 
10.2 < log{M) < 10.6 
10.6 < log{M) < 11.0 
11.0 < log{M) < 11.4 
11.4 < log{M) < 11.8 


-3.81 ±0.28 
-3.86±0.18 
-3.81 ±0.10 
-3.83 ±0.09 
-4.13 ±0.09 
-4.74 ±0.18 


-3.51 ±0.09 
-3.57 ±0.11 
-3.25 ±0.10 
-3.45 ±0.12 
-3.79 ±0.17 


-3.61 ±0.08 
-3.59 ±0.09 
-3.77 ±0.11 


-3.34 ±0.08 
-3.31 ±0.08 
-3.95 ±0.14 


-3.24 ±0.07 
-3.55 ±0.10 
-4.75 ±0.21 



in close major-merger pairs in this work and in Bundy et al. 
(2009). The Y parameter derived using paired galaxies in 
our sample (the inverse of Eq. 11) and that of Conselice et al. 
(2009) derived using morphologically selected mergers are in 
very good agreement. The higher T values of Bundy et al. 
(2009) are likely due to the relatively long merging time scale 
they adopted from Kitzbichler & White (2008). On the other 
hand, the low T values and lack of evolution of Lotz et al. 
(2008) are because of the inclusion of minor mergers in their 
sample. Our results are in good agreement with the predic- 
tion of the default (semi-empirical) model of Hopkins et al. 
(2010b). 

6. MAJOR MERGERS, ELLIPTICAL GALAXY FORMATION, AND 
GALAXY ASSEMBLY 

6. 1 . Mass Dependent Volume Merger Rate 

The volume merger rate, Ry, measures the frequency of 
merger events in a given volume in the universe. Here we 
define the mass dependent Ry in terms of the stellar mass of 
the merger remnant, which is the total stellar mass of the two 
galaxies involved in the merging (ignoring the mass of stars 
formed during the merger): 

Rv(M star ,z) = 0.5 x R mg (M star /lO o - 2 ,z)0(M star /lO a2 ,z), 

(12) 

where 0(M,z) is the GSMF of galaxies in the parent sam- 
ple (D09), the factor of 0.5 is due to the fact that every ma- 
jor merger event involves two galaxies of similar mass. We 
also assume that on average the mass of a merger remnant is 
0.2 dex higher than that of individual galaxies involved in the 
merger. This is because, under the assumption that mass ra- 
tio distribution is flat (Appendix C), pairs in our sample has 
a mean mass ratio of 0.2 dex. Therefore the logarithm of the 
mean ratio between the total mass of a pair and that of the 
primary is log(l + 10~° 2 ) ~ 0.2. Our results on the Ry are 
presented in Table 6. 

In Fig. 8 we compare our results with the average growth 
rate of elliptical galaxies (dashed line) and that of RQGs (dot- 
ted line) between z=0.3 - 0.9, estimated using the differ- 
ences between their GSMFs at z=0.3 and z=0.9, taken from 
Ilbert et al. (2010), divided by 3.88 Gyr (the time span corre- 
sponding to the redshift interval of [0.3,0.9]). The solid line 
is the average volume merger rate calculated using Eq. 12 by 
replacing R mg with its best fit (Eq. 11) and averaged over 
the same redshift range of 0.3 - 0.9. It shows that major 
mergers can fully account for the formation of both massive 
ellipticals and RQGs (M star > 10 10 5 M ). This contradicts 
Bundy et al. (2009) who concluded that the major-merger rate 
is too low to fully explain the formation of RQGs since z=l. 
The major reason for the contradiction is due to the differ- 
ence in the adopted merger time scales in this work and in 



Bundy et al. (2009): Our T mg , derived from the results of 
Lotz et al. (2010), is about a factor of 2 shorter than that used 
by Bundy et al. (2009). There is also a difference in the for- 
mation rates of ellipticals and RQGs adopted in this work (es- 
timated from results of Ilbert et al. (2010)) and in Bundy et al. 
(2009). The latter is about 50 - 100% higher than the former. 

We define "dry mergers" ("wet mergers") as those in pairs 
or multiple systems consisted of only RQGs (SFGs), and 
"mixed mergers" the rest of galaxies in the pair sample. In 
Fig. 9 we compare the volume merger rates of dry merg- 
ers, and those of wet and mixed mergers combined, to the 
formation of ellipticals. Our results show that wet/mixed 
mergers alone can account for the formation rate of mas- 
sive ellipticals and RQGs, even for the most massive ones of 
M s tar > 10 11 - 3 M . Our results are consistent with Lin et al. 
(2008), who also found that the wet and mixed mergers dom- 
inated over the dry mergers since z~ 1 . 

It should be pointed out that ellipticals (Es) and 
RQGs are not identical (albeit with large overlap) since 
there are both red disk galaxies (Bamford et al. 2009; 
Bundy et al. 2010) and blue ellipticals (Kannappan et al. 
2009; Huertas-Company et al. 2010). Both Es and RQGs 
are likely originated from star forming disk galaxies: RQGs 
formed through SFR quenching (Bell et al. 2007; Faber et al. 
2007), and ellipticals through mergers (Toomre 1978; Barnes 
1988) or secular evolution (Kormendy & Kennicutt 2004). 
Hopkins et al. (2008) argue that only major mergers can 
reproduce the kinematic properties of massive ellipticals, 
whereas some low mass ellipticals are "pseudo-bulges" 
formed through disk instabilities and secular evolution in late 
type galaxies. 

Many quenching mechanisms in the literature are related 
either directly to merger induced feed-backs (e.g. gas con- 
sumption by extreme starbursts and gas loss due to super- 
winds) or to massive bulges (such as the AGN quench- 
ing, Bensen et al. 2003; Somerville et al. 2008, and morpho- 
logical quenching, Martig et al. 2009). Hence a close re- 
lation between major mergers and RQGs formation is ex- 
pected. Indeed Hopkins et al. (2008) argued that a wide range 
of observations (e.g. the bivariate red fraction as a func- 
tion of galaxy and halo mass, the density of passive galax- 
ies at high redshifts, the emergence/evolution of the color- 
morphology-density relations at high redshift, and the fraction 
of disky/boxy spheroids as a function of mass) favor a major- 
merger related quenching model to other quenching models. 

Our results (Fig. 8 and Fig. 9) support the model of 
Hopkins et al. (2008). Major mergers, dominantly wet or 
mixed, can fully account for the formation rates of both el- 
lipticals and RQGs with M sta r > 10 10 5 M . For most mas- 
sive galaxies with M star > 10 1L3 M , the major-merger rate 
agrees very well with the two formation rates (in this mass 
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Figure 9. Left: Mass-dependent volume merger rates of wet (S+S) and mixed (S+E) mergers, in different photo-z bins. The dashed line is the 
average growth rate of ellipticals between z=0.3 - 0.9 and dotted line that of RQGs, both being estimated using results of Ilbert et al. (2010). 
Right: Mass-dependent volume merger rates of dry (E+E) mergers. 
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range most ellipticals and RQGs belong to the same popula- 
tion of red elliptical galaxies). In the mass range of 10 10 5 < 
M st ar < 10 11 - 3 M , the major-merger rate is slightly higher 
than both formation rates. Two factors may be responsible for 
this: (1) Remnants of some gas-rich wet mergers may remain 
to be blue disk galaxies (Hopkins et al. 2009). (2) Dry merg- 
ers, contributing most in this mass range, may move some red 
ellipticals to higher mass. 

Fig. 8 and Fig. 9 also show that most (i.e. ~ 2/3) of 
low mass ellipticals and RQGs (M star £ 10 10 3 M ) are not 

produced by mergers. Many authors have argued that these 
galaxies are mostly quenched satellite galaxies whose gas 
halos are stripped by much more massive central galaxies 
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Figure 11. The volume rates of (U)LIRGs (10 1L5 < Li R /L < 10 12 5 ), in 
different redshift and stellar mass bins. Estimated using data taken from 
Kartaltepe et al. (2010) and assuming a (U)LIRG time scale of 140 Myr. The 
solid line is the average mass dependent volume merger rate (this work), iden- 
tical to that in Fig. 8. 



(van den Bosch et al. 2008; Peng et al. 2010). 

6.2. Impacts of Major Mergers on Galaxy Assembly 

Mergers shift galaxies from lower mass bins to higher mass 
bins in the GSMF. The efficiency of this process is clearly 
critical for the hierarchical structure formation paradigm. 
Drory & Alvares (2008) tried to answer this question via com- 
parisons between observed GSMF variation against redshift 
and that predicted by the SFR vs. z relation which was well 
established (in particular for z < 1), attributing the difference 
to the merger effects (including both major and minor merg- 
ers). Our results on major-merger rates provide a new and 
more direct approach. 

to evaluate the 



We define the following parameter 



'merger 



Galaxy Pairs in COSMOS - Merger Rate Evolution Since z = 1 



11 



impact of major mergers on the GSMF 

0merger(M star ) = Rv(M star )/0(M star ) - R mg (M star ), (13) 

where 0(M star ) is the GSMF (D09), R v (M star ) the volume 
merger rate defined in Eq. 12, and R mg the differential merger 
rate estimated using Eq. 11. 

In Fig. 10 we compare our results with those of 
Drory & Alvares (2008). The solid line is our result on the 
mean (^merger (M star ) over 0.3 < z < 0.9. In this redshift range, 
major mergers (as opposed to star formation) have significant 
impact to the galaxy mass assembly only for the most mas- 
sive galaxies with M star ^ 10 1L3 . The GSMF change due to 
major mergers dominates that due to star formation only at 
M star ^ 10 11 - 3 . For less massive galaxies with M star < 10 112 
the GSMF change due to major mergers is negligible (ampli- 
tude < 10 % Gyr -1 ) in comparison to that due to star for- 
mation. For these galaxies, our result is much flatter than 
that of Drory & Alvares (2008) for the GSMF change rate 
due to mergers at both z=0.5 and z=l. For massive galax- 
ies (M s tar ~ 10 1L3 ) we find much steeper mass dependence 

than Drory & Alvares (2008). The major reason of the dis- 
crepancy is due to the difference between the redshift depen- 
dent GMSFs used in this work (Drory et al. 2009) and those 
in Drory & Alvares (2008), the latter were derived using data 
from earlier FDF/GOODS surveys. 

7. MAJOR MERGERS AND (U)LIRGS 

Kartaltepe et al. (2010) carried out a study of luminous 
IR galaxies (LIRGs: 10 11 < L IR /L < 10 12 ) and ultra- 
luminous IR galaxies (ULIRGs: L lR /L > 10 12 ) in the S- 
COSMOS survey (Sanders et al. 2007). We neglect ULIRGs 
of L IR /L > 10 12 - 5 since very few galaxies have such high 
IR luminosities, and the AGN fraction increases rapidly with 
the L IR among these galaxies (Kartaltepe et al. 2010). Taking 
galaxies with 10 1L5 < L lR /L < 10 12 5 from their sample and 
adopting a (U)LIRG time scale of 140 Myr (Kartaltepe et al. 
2010), we estimate the (U)LIRG rates in different stellar mass 
bins and in the redshift range of 0.2 < z < 1.0. The results are 
presented in Fig. 11, compared with the average mass depen- 
dent volume merger rate. 

The mass dependence of (U)LIRG rates in all photo-z bins 
have the shape of the log-normal function, peaking at a rather 
constant mass of ~ log(M* tar /M®) =10.8. In the low photo- 
z bin (0.2 < z < 0.4) (U)LIRGs are less frequent, consis- 
tent with the fact that (U)LIRGs are very rare in the lo- 
cal universe (Sanders & Mirabel 1996). The average mass 
dependent volume merger rate is above or comparable to 
the (U)LIRG rates in all redshift and mass bins, and there- 
fore it is possible that most of the (U)LIRGs in this red- 
shift range are major-mergers, just like their local counter- 
parts (Sanders & Mirabel 1996). Using morphological classi- 
fications, Kartaltepe et al. (2010) found that £ 50% of their 

(U)LIRGs are major-mergers. This means that the merger- 
induced (U)LIRG rates are even more below the average 
merger rate, in particular for M sta r ~ 10 10 3 andM sta r ~ 10 1L3 , 

than being depicted in Fig. 11. Hence, it is likely that 
a large fraction of major-mergers, in particular those with 
Mstar £ 10 10 - 3 or Mstar £ 10 1L3 , may not become (U)LIRGs. 
It is interesting to note that the most massive mergers of 
Mstar £ 10 1L3 M have rather low (U)LIRG rate. Most of 



them are wet or mixed mergers, but many probably have rel- 
atively low gas content. Galaxies of lower mass (M sta r °° 
10 10 M ) also have low (U)LIRG rate because their gas mass 
is not adequate to sustain the very high SFR of the extreme 
starbursts in (U)LIRGs. 

8. SUMMARY 

We have presented results of a statistical study on the cos- 
mic evolution of the mass dependent major-merger rate since 
z = 1. A stellar mass limited sample of major-merger pairs 
(the CPAIR sample) was selected from the archive of the 
COSMOS survey. It includes 617 galaxies in pairs/multiple- 
systems with stellar mass ratios /ileq2.5, projected separa- 
tions in the range of 5 < r pro j < 20 fr 1 kpc, and in the photo-z 
range of 0.2 < z pri0 t < 1.0. The pair selection was based on 
photo-z, with the criterion of Az p hot / (1 + z P hot) < 0.03, and 
on visual inspections of the HST-ACS images. The CPAIR 
sample is divided into four photo-z bins of [0.2 < z < 0.4, 
0.4 < z < 0.6, 0.6 < z < 0.8, 0.8 < z < 1]. Various biases 
in the sample selection that caused incompleteness and spu- 
rious pair contaminations have been studied. This resulted in 
a completeness correction factor of [0.78 ± 0.05, 0.74 ± 0.05, 
0.71 ± 0.05, 0.60 ± 0.05] and a reliability correction factor of 
[0.80 ± 0.06, 0.79 ± 0.06, 0.77 ± 0.06, 0.78 ± 0.06], respec- 
tively, for pairs in the four photo-z bins. The CPAIR sample 
is complemented by a local (z=0) major-merger pair sample, 
selected in the K-band from cross matches between 2MASS 
and SDSS-DR5 (an updated version of the KPAIR sample in 
Domingue et al. 2009). 

Mass dependent pair fractions at different redshifts derived 
using CPAIR and KPAIR samples show no significant varia- 
tions with stellar mass. The integral pair fraction (i.e. pair 
fraction of all galaxies regardless of stellar mass) demon- 
strates a moderately strong cosmic evolution, with the best- 
fitting function of fpair = 10- 1 - 88 ( ±0 - 03 )(l+z) 2 - 2 ( ±0 - 2 ). 

The merger time scale was taken from the simu- 
lation results of Lotz et al. (2010): T mg /Gyr = 0.3 x 
(M s tar/lO lo - 7 M )- a3 (l +z/8). The best-fitting function 
for the differential merger rate (for ji < 3 mergers) is 
Rmg/Gyr 1 =0.0453 x (M s tar/lO ia7 M )°- 3 (l+z) 2 - 2 /(l+z/8). 
Accordingly, on average, galaxies of M sta r ~ 1 10-1 1-5 M 
have undergone ~ 0.5 - 1.5 times major mergers since z=l. 
Our result on the differential merger rate is in very good 
agreement with that estimated using morphologically selected 
major-mergers (Conselice et al. 2009) and with the prediction 
of the "semi-empirical" model of Hopkins et al. (2010b). 

The mass dependent major-merger rates derived in this 
work indicate that, for massive galaxies (M sta r > 10 10 5 M ) 
at z < 1, major mergers involving star forming galaxies (i.e. 
wet and mixed mergers) can fully account for the formation 
rates of both ellipticals and red quiescent galaxies (RQGs), 
lending support to models that link both bulge formation and 
SFR quenching to major mergers (e.g. Hopkins et al. 2008). 
On the other hand, major mergers cannot be responsible for 
the formation of most low mass ellipticals and RQGs of 
Mstar £ 10 10 - 3 M . Dry mergers contribute negligibly to the 

major-merger rate in all mass and photo-z bins. Major merg- 
ers have significant impact to the stellar mass assembly of the 
most massive galaxies (M sta r > 10 1L3 M ). For less massive 
galaxies the stellar mass assembly is dominated by the star 
formation. 

Comparisons with mass dependent (U)LIRG rates in differ- 
ent redshift bins suggest that the frequency of major-merger 
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events is comparable or higher than that of (U)LIRGs. Most 
low mass mergers (M star < 10 10 3 M ) and most very massive 

mergers (M star £ 1O 1L3 M ) may not become (U)LIRGs. 
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APPENDIX 

A. INCOMPLETENESS DUE TO MISSING VERY CLOSE PAIRS — ANALYSIS 
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Figure A-l. Plot of the angular separation e vs. redshift relation for pair candidates (open diamonds), and for very close ACS pairs that are missed by the photo-z 
pair candidate sample (plus symbols). 



In order to estimate how many pairs with e ^ 2" are missing in our sample, we carried out an analysis exploiting the COSMOS 

HST-ACS lensing catalog (Leauthaud et al. 2007, 2010). It includes 1.2 x 10 6 galaxies detected by HST-ACS in the F814 band 
(hereafter i A cs band), with an angular resolution of FWHM = 0.12" (Leauthaud et al. 2007). From this catalog, we selected a 
sample of "very close ACS pairs" through the following procedure: 

: (1) Find the match in the HST-ACS lensing catalog for every D09 galaxy of log(M star ) > log(Mii m ) (see Section 2 for the 
definition of Mii m ) with a searching radius of 0.8" and the criterion of |/acs -*ab,do9| < 1 mag. 

: (2) Around the ACS matches of D09 galaxies, we search for ACS pairs with three criteria: (i) | Az'acsI < 1 mag; (ii) e < 2.0"; 
(iii) r pr0 j > 5 IT 1 kpc. 
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Figure A-2. Comparison of mass distributions of CPAIR galaxies and of galaxies in the single photo-z ACS pairs (e < 2") with 0.8 < z p h t < 1-0, including 31 
of the 47 "very close" ACS pairs. 

The procedure selected 222 very close ACS pairs. Among them, 171 pairs have both component galaxies with photo-z matches, 
including 53 pairs found in the photo-z pair sample and the remaining 118 pairs consisting of galaxies of discordant photo-z 's or 
with mass ratios larger than 2.5. In the remaining 51 pairs, 4 are multi-peak single galaxies in the ACS images. The final sample 
has 47 pairs in which only one of two galaxies was detected in the photo-z catalog. These pairs, shown in Fig. A-l by red crosses, 
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have the average angular separation < e >= 1.45" with the standard deviation of a = 0.21" . Fig. A-l also shows that, for pairs of 
z > 0.5, the lower boundary for the pair separations, r pro j > 5 fr 1 kpc, corresponds to a angular separation of e ~ \" . 

We then used Monte Carlo simulations to estimate the expected number of spurious pairs in the sample of very close ACS 
pairs, utilizing the 138001 galaxies in D09 sample. In each of the 100 simulations, we randomly put these galaxies in a 1.7 deg 2 
region, with all other properties of the galaxies, including the photo-z and stellar mass, intact. We then search companions around 
each of the galaxies of log(M sta r) > log(Mii m ) in the simulated sample according to the following criteria: |A/ab| < 1 mag, (ii) 
e < 2.0", and (iii) r pro j > 5 h _1 kpc. Spurious pairs that do not pass the four pair selection criteria in Section 2 were then counted. 
These simulations found a mean spurious pair number of 103.3 with a l-a dispersion of 7.0. As described above, the number of 
confirmed spurious pairs of photo-z galaxies in the sample of very close ACS pairs is 118. This number is slightly higher than 
the mean total number of spurious pairs (103.3±7.0) predicted by the Monte Carlo simulations, perhaps due to galaxy clustering. 
Thus majority of 47 single photo-z ACS pairs are in fact real and not chance superpositions. Indeed their ACS images very often 
show signs of interaction. In what follows we shall make the conservative assumption that all 47 ACS pairs are real major-merger 
pairs. 

There are Nacs = [1,5, 10,31] single photo-z ACS pairs in the four photo-z bins. The corresponding incompleteness due to 
missing of very close pairs, estimated according to the ratio Nacs /(N p h to-z + Nacs) (N p h to-z being the number of photo-z pair 
candidates), is 0.01 ±0.01, 0.06±0.03, 0.08±0.03, and 0.20±0.04 for the four photo-z bins, respectively. As a check, we found 
no significant difference between the stellar mass distributions of galaxies in the ACS single photo-z pairs and of those in the 
CPAIR sample in the redshift bin of 0.8 < z phot < 1 .0 (Fig. A-2). 

B. INCOMPLETENESS DUE TO PHOTO-Z ERRORS AND SPURIOUS PAIRS DUE TO PROJECTION — MONTE CARLO SIMULATIONS 
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We examine the accuracy of the photo-z using the spectroscopically measured redshift (hereafter spec-z) of galaxies in the 
parent sample that were observed in the zCOSMOS survey (Lilly et al. 2007). zCOSMOS includes a magnitude-limited survey 
(zCOSMOS-bright) for about 20,000 galaxies of / AB < 22.5 and 0.1 < z < 1.2, covering 1.7 deg 2 COSMOS includes 10643 
galaxies. There are 5001 matches (matching radius = 2") between galaxies in the parent sample (Table 1) and zCOSMOS sources 
with reliable spec-z measurements (z-class indices being 4's, 3's, 9.5, 9.4, 9.3, 2.5, 2.4 or 1.5). This is 19.5% of 2571 1 galaxies in 
the parent sample that have /ab < 22.5. Fig. B-l shows (z p hot-z spec )/(l +z spec ) of the 5001 galaxies. The distribution, an estimate 
of the photo-z error distribution (the spec-z error ;S 100 km sec -1 , Lilly et al. 2007), has a 1 .48 x median(|z p hot-z spec | / (l +z spec )) = 
0.007, consistent with the result of Ilbert et al. (2009). It can be best fitted by a Lorentzian function of the form 



L(x) = 



7?r l + ((x-x )/7) 2 ' 



(Bl) 



with x = (z phot - z spec ) /(l + z spec ), 7 = 0.0056 ± 0.0004 and x = 0.0005 ± 0.0003. 

In Fig. B-2, the i-band magnitudes are shown for galaxies in pair candidates (red squares) and in the parent sample (black dots). 
Both the pair sample and the parent sample are dominated by galaxies brighter than /ab = 22.5. In the pair candidates, factions of 
0.032, 0.061, 0.168, 0.352 of the sample are fainter than / AB = 22.5 for the four redshift bins. 

The Monte Carlo simulations for the estimation of completeness correction factor due to photo-z errors were based on the 
above results. The simulations included a total of 1000 repeats for each photo-z bin, each containing 100 pairs. For pairs in a 
given photo-z bin, each galaxy had a chance to be fainter (or brighter) than /ab = 22.5, the probability (simulated by a random 
number generator) being equal to the observed fraction of such galaxies in the photo-z bin. For galaxies brighter than /ab = 22.5, 
an error of z p h t/(l + z p hot) was assigned to it using a random number generator with weighted probability distribution function 
given by a Lorentzian with xo = and 7 = 0.0056. For galaxies fainter than /ab = 22.5, the random error was generated according 
to a Lorentzian function with x = and 7 = 0.0096 (~ 0.0056 x 0.012/0.007, where 0.012/0.007 is the ratio between the photo- 
z accuracies for galaxies of 22.5 < /ab < 24 and of /ab < 22.5, Ilbert et al. 2009). We ignored the real velocity difference and 
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Figure B-2. The i-band magnitude (z'ab) vs. photo-z plot. Galaxies in the pair candidates are shown by red squares, and galaxies in the parent sample (Table 1) 
by black dots. The dashed line marks the boundary of z'ab =22.5. 

assumed that the photo-z difference in a pair is purely due to photo-z errors. The completeness factor was estimated by the 
fraction of simulated pairs with | Az p h t/(1 +z p hot)| < 0.03. According to the simulations, the completeness correction factor is 
[0.78 ± 0.05, 0.77 ± 0.06, 0.76 ± 0.05, 0.73 ± 0.05] in the four photo-z bins. 

Another set of Monte Carlo simulations (each consisting of 1000 repeats) were carried out to estimate the SPF due to projection, 
utilizing the 138001 galaxies in D09 sample. Here we pretended that the photo-z's in that sample are 100% accurate, and then 
added errors to them using the same algorithm as described above. The sky coordinates of the galaxies were also randomized, 
filling a 1.7 deg 2 region uniformally. The other properties of the galaxies, including the stellar mass, were left intact. We then 
selected pairs from this simulated parent sample by applying the four selection criteria presented in Section 2. Spurious pairs 
were identified when the "true velocity difference", calculated using the "true redshifts" (i.e. the photo-z's without added error), 
is > 500 km sec -1 . In the four photo-z bins, the simulations found [7.5 ± 2.8, 6.2 ± 2.5, 11.4 ± 3.3, 11.1 ± 3.4] spurious pairs. 
Dividing these numbers by the numbers of photo-z pair candidates, the predicted SPF in the four photo-z bins are 0.07 ± 0.03, 
0.08 ± 0.03, 0. 10 ± 0.03 and 0.09 ± 0.03, respectively. 

C. DISTRIBUTION OF LINE-OF-SIGHT VELOCITY DIFFERENCE OF MAJOR-MERGER PAIRS 
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Figure C-l. Cumulative distribution of Av of major-merger pairs, derived using a sample of z=0 K-band selected pairs (Domingue et al. 2009). 
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Local major-merger pairs in Domingue et al. (2009) were selected using nearly identical selection criteria as the CPAIR sample 
(cf. Section 4) except for that spec-z were used and they include all pairs (isolated or in groups/clusters) with Av spec < 1000 
km sec -1 . Here we exploit these data to determine the cumulative distribution of Av and, in particular, the fraction of pairs with 
500 < Av < 1000 km sec -1 . We assume that all pairs with Av > 1000 km sec -1 are spurious. In the sample of Domingue et al. 
(2009), which is 0.82% complete for the spec-z, 135 pairs have measured v spe c for both components and v spe c > 2000 km sec -1 
(excluding pairs in the local super-cluster). The cumulative distribution of Av of these pairs is plotted in Fig. C-l. From this 
distribution, we found that the fraction of pairs with Av > 500 km sec -1 is 9.6% with a random error of 2.5% (binomial statistics). 
It should be pointed out that the Av distribution is sensitive to the environment and to the pair separation (Ellison et al. 2010). 
Therefore caution should be taken when applying the result here to other pair samples. 

D. INCOMPLETENESS AND SPURIOUS PAIRS FRACTION — A COMPARISON WITH SPEC-Z PAIRS 

We made a comparison between pairs in the CPAIR sample and a sample of spec-z pairs (ZPAIR sample) selected from 
the zCOSMOS survey (Lilly et al. 2007). In principle, this comparison applies only to galaxies brighter than z A b = 22.5, the 
magnitude limit of zCOSMOS. However, most of galaxies in photo-z pairs are brighter than /ab = 22.5 (Fig. B-2). Images of 
the 16 pairs in the ZPAIR sample, listed in Table D-l, are shown in Fig. D-l. Among them are all CPAIRs (14) where both 
components have spec-z's. They comprise a very small fraction of total CPAIR sample because only 19.5% of galaxies of 
/ab < 22.5 in the parent sample have spec-z. All but two of the 14 pairs have the velocity difference Av spe c < 500 km sec -1 . This 
results in an empirical estimate for the SPF of 2/14 = 0.14, or a reliability of 0.86, with a binomial uncertainty of ±0.12. 

Table D-l also includes two pairs (ZPAIR-03 and ZPAIR-07) that were missed by CPAIR sample, both having Az p h t/(1 + 
Zphot) > 0-03 but Av spec < 500 km sec -1 . As shown in Fig. D-l, the primary of ZPAIR-03 is a close pair itself, and its large z p h t 
error (Table D-l) is likely due to confusion in the photometric data. ZPAIR-07, having Az p h t/(1 +z p hot) = 0.031, barely missed 
the cut for the photo-z pairs. Adding the 12 genuine CPAIRs, there are 14 pairs with Av < 500 km sec -1 in the ZPAIR sample. 
This results in a completeness for the CPAIR sample of 12/14 = 0.86, with a statistical error 0.12. 

E. DISTRIBUTION OF PRIMARY-TO-SECONDARY MASS RATIO OF CLOSE PAIRS 
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Figure E-l. Left: Differential distributions of log(/j) (/j, = M pr i/M2 n d) of close pairs (5 < r pro j < 20 h 1 kpc) in three mass bins. Derived using pair candidates 
selected in a volume limited sample of galaxies with redshift in the range of 0.2 < z < 0.4. Right: Normalized cumulative distributions of log(/i) (/j = M pr i/M2 n d) 
of the same close pairs. 



In the parent sample of CPAIR, galaxies in the photo-z bin of 0.2 < z p h t < 0.4 form a volume-limited sample above the 
stellar mass limit M star = 10 9 °M . Using this volume limited sample and applying the same pair selection criteria in Section 2 
except for expanding the mass ratio limit to log(/i max ) = 1.0, where ji = M pri /M2 n d is the mass ratio, we selected pair candidates 
(including both major and minor pairs) in three mass bins: 10.0 < log(Mp r i/M ) < 10.4, 10.4 < log(M pr i/M ) < 10.8, and 
10.8 < log(Mp r i/M ) < 1 1.2. Assuming that completeness and reliability corrections for these pair candidates do not depend on 
the mass ratio /i, we calculated the differential and cumulative log(/i) distributions of these close pairs (5 < r pro j < 20 fT 1 kpc). The 
results are plotted in Fig. E-l. It shows that the flat distribution, i.e. df pa i r /dlog(/i) = constant, is a reasonably good approximation. 
It is worth noting that our result is different from that of Ellison et al. (2010), who found a mass ratio distribution for SDSS pairs 
that is tilted toward low fi pairs (i.e. major mergers). However, their result is affected significantly by the "missing secondary" 
bias, causing severe incompleteness of the ji ~ 10 minor-mergers in their sample. 




Figure D-l. HST-ACS (F814) images of close major-merger pairs in zCOSMOS survey (ZPAIRs in Table D-l). The size of all images is 20" x 20". The 
crosses mark the positions of component galaxies in the pairs. Notes: ZPAIR-03 and ZPAIR-07 are not in the CPAIR (photo-z pairs) sample. The large neighbor 
in ZPAIR-14 is likely a foreground galaxy. ZPAIR-16 is in a compact group of galaxies (three in this group were included in the CPAIR sample). 
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Table D-l 

Pairs in Z-COSMOS Survey (ZPAIR) 



ZPAIR 


Av S pec 


RAi 


Deci 


z 1 


zclassi 


Zphot,l 


l0g(M sta r,l) 


ID 


(km sec -1 ) 


(degree) 


(degree) 








(M )' 




t—i^phot/ \ L ~ ^pnoU 


RA 2 


Dec2 


Zspec,2 


zclass2 


Zphot,2 


lw &V 1VA star,2/ 






(degree) 


(degree) 








(M©) 


01 


43 


149.713610 


2.019610 


0.6444 


4.5 


0.6518 


10.92 




0.012 


149.713950 


2.019773 


0.6441 


4.5 


0.6324 


10.83 


02 


218 


149.839700 


1.929228 


0.3722 


3.5 


0.3856 


10.80 




0.020 


149.838780 


1.930125 


0.3711 


4.5 


0.3585 


10.76 


03 t 


70 


150.009050 


2.274964 


0.4726 


2.5 


0.5769 


11.05 




0.071 


150.008210 


2.275954 


0.4730 


2.5 


0.4645 


10.78 


04* 


817 


150.107580 


2.556516 


0.5038 


3.5 


0.4915 


10.33 




0.004 


150.107510 


2.557509 


0.4990 


3.5 


0.4969 


10.21 


05 


73 


150.115480 


1.975120 


0.4385 


4.5 


0.4431 


10.45 




0.025 


150.115680 


1.976216 


0.4381 


3.5 


0.4069 


10.05 


06 


65 


150.126060 


1.913758 


0.7360 


2.5 


0.7168 


11.17 




0.003 


150.125690 


1.913726 


0.7365 


3.5 


0.7108 


11.15 


07 t 


11 


150.168760 


2.315481 


0.8524 


2.5 


0.7921 


10.97 




0.031 


150.168890 


2.316234 


0.8523 


1.5 


0.8473 


10.94 


08 


220 


150.196380 


2.371582 


0.6834 


4.5 


0.6783 


10.86 




0.003 


150.196460 


2.370591 


0.6850 


22.5 


0.6726 


10.64 


09 


103 


150.230880 


1.845002 


0.6226 


2.5 


0.6072 


10.67 




0.014 


150.230550 


1.844713 


0.6233 


3.5 


0.5840 


10.60 


10 


118 


150.258800 


1.988773 


0.7258 


2.5 


0.7168 


10.70 




0.001 


150.258510 


1.988547 


0.7267 


2.5 


0.7191 


10.32 


11* 


6318 


150.359560 


2.659517 


0.4309 


22.5 


0.3991 


10.94 




0.002 


150.359280 


2.660543 


0.3974 


3.5 


0.4018 


10.81 


12 


142 


150.396810 


2.519130 


0.2189 


24.5 


0.2283 


9.97 




0.003 


150.396980 


2.517332 


0.2195 


4.5 


0.2323 


9.93 


13 


333 


150.421160 


2.654335 


0.2144 


4.5 


0.2312 


9.57 




0.003 


150.421480 


2.654870 


0.2158 


4.5 


0.2271 


9.21 


14 


118 


150.457210 


2.695287 


0.2189 


4.5 


0.2277 


9.45 




0.008 


150.456970 


2.696160 


0.2194 


4.5 


0.2184 


9.07 


15 


219 


150.494640 


2.187936 


0.3695 


4.5 


0.3674 


9.97 




0.007 


150.494120 


2.187388 


0.3684 


3.5 


0.3772 


9.58 


16 


196 


150.504990 


2.225083 


0.8374 


22.5 


0.8360 


11.72 




0.002 


150.504730 


2.224246 


0.8357 


2.5 


0.8315 


11.34 



Note. 

f Missing in the CPAIR sample. 

* Spurious pairs (with Av spec > 500 km sec~ l ) in the CPAIR sample. 



